Symmetric Versus Nonsymmetric Structure of the Phosphorus Vacancy on InP(llO) 
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The atomic and electronic structure of positively charged P vacancies on InP(llO) surfaces is 
determined by combining scanning tunneling microscopy, photoelectron spectroscopy, and density- 
functional theory calculations. The vacancy exhibits a nonsymmetric rebonded atomic configuration 
with a charge transfer level 0.75 ±0.1 eV above the valence band maximum. The scanning tunneling 
microscopy (STM) images show only a time average of two degenerate geometries, due to a thermal 
flip motion between the mirror configurations. This leads to an apparently symmetric STM image, 
although the ground state atomic structure is nonsymmetric. Copyright 2000 by the American 
Physical Society. 
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Although it is well known that point defects can exert a 
profound influence on the physical properties of semicon- 
ductors, the determination of the atomic scale geometric 
and electronic structure of point defects has remained a 
challenging task for theoretical as well as experimental re- 
search. One particularly striking example is anion vacan- 
cies on (110) surfaces of III-V semiconductors, where no 
agreement has been reached regarding even basic proper- 
ties, such as (i) the symmetry of the atomic structure and 
(ii) the energy of defect levels in the band gap: Although 
recent density-functional theory (DFT) calculations for 
positively charged arsenic vacancies on the GaAs (110) 
surface agreed that the gallium atoms neighboring the va- 
cancy relax into the surface layer (in contrast with an ear- 
lier tight-binding calculation (!]]), one calculation found a 
rebonded configuration breaking the mirror symmetry of 
the surface ||, whereas the other predicted a fully sym- 
metric configuration to have the lowest energy ||. High 
resolution scanning tunneling microscopy (STM) images 
show a density of states preserving the mirror symme- 
try of the surface at the defect site |Q,H|- Although this 
seems to favor a symmetric atomic structure, the exper- 
imental results could not be matched to results of any 
of the DFT calculations |^-||. Furthermore, scanning 
tunneling spectroscopy (STS) yielded a local downward 
band bending of 0.1 eV [Q, whereas surface photovolt- 
age measurements found a band bending of 0.53 ± 0.3 eV 
H at the site of the positively charged As vacancy on p- 
doped GaAs(llO) surfaces. Concerning the energy levels, 
the two different DFT calculations predicted the charge 
transfer levels (+/0) to be 0.32 eV g and 0.1 eV @, 
and the lowest Kohn-Sham eigenvalues in the band gap 
to be 0.73 eV @ and 0.06 eV above the valence band 
maximum (VBM) . In view of this puzzling situation it is 
obvious that the interpretation of the experimental STM 
images as well as the structure of the anion vacancies 
on (110) surfaces of III-V semiconductors are still under 
debate 



Discrepancies such as those pointed out above can arise 
due to limitations of the methods used. On the theoreti- 
cal side DFT calculations have proven to be powerful to 
determine the structure of point defects [|l0|. Yet differ- 
ences in the size of the supercell, the pseudopotentials 
employed, or the exchange-correlation functional imple- 
mented might lead to different results, such as the de- 
fect's symmetry or the position of the defect levels, espe- 
cially in situations where two configurations are almost 
degenerate in energy. On the experimental side the STM 
images alone do not provide direct information about the 
atomic structure nor are the experimental conditions, no- 
tably the concentration of the vacancies, the same in all 
experiments. Furthermore the quantities discussed to de- 
scribe the position of the defect level, i.e., Kohn-Sham 
eigenvalues within DFT, band bending, and charge trans- 
fer levels, are different physical quantities and thus (typ- 
ically) deviate from each other significantly. In addition, 
the tip of the STM may affect the structure of the va- 
cancy, since it is well known that the tip can even excite 
vacancies to migrate pj . Therefore the combination of 
state-of-the-art theory and well defined experiments is 
needed to get reliable results and test the correct inter- 
pretation of each of the separate results. 

In this letter we combine three different methods 
to determine the geometric and electronic structure of 
phosphorus (P) vacancies on p-doped InP(llO) surfaces: 
STM, photoelectron spectroscopy (PES), and DFT. This 
enables us to determine the charge transfer level of the 
defect state in the band gap and the equilibrium atomic 
structure of the vacancy in p-type material. We identify 
the rebonded structure as the ground state of the pos- 
itively charged defect and demonstrate that a thermal 
flip motion between two mirror configurations results in 
simulated STM images consistent with the experimen- 
tal findings. Our results solve the apparent discrepancies 
of the previous experimental and theoretical studies [O , 
provide a comprehensive picture of the vacancy structure, 
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FIG. 1. Occupied state STM image of P monovacancies 
formed by annealing p-doped InP(llO) surfaces at 160° C. The 
image has been acquired at —2.7 V tunneling voltage and 0.9 
nA current. 



and point out a methodology to determine the structure 
of defects with high accuracy. 

The surface vacancies were produced in high con- 
centrations by annealing initially defect-free p-doped 
InP(llO) cleavage surfaces [carrier concentration n^op of 
(1.3-2.1) x 10 18 cm- 3 ] at 160°C @Q. Figure 1 shows 
that this procedure leads indeed to only one type of de- 
fect on the surface which was identified as a singly posi- 
tively charged phosphorus vacancy previously and whose 
concentration can be controlled by the annealing time 
JLsl p^[ . We prepared several samples from the same 
wafer with a variety of annealing times in an ultrahigh 
vacuum (UHV) chamber at the " Berliner Elektronenspe- 
icherring fur Synchrotronstrahlung" storage ring for the 
photoelectron spectroscopy and in another UHV system 
for the STM measurements ensuring an exact control of 
defect concentration in both experiments. All measure- 
ments were done at room temperature. Angle-resolved 
PES measurements were used to probe the valence band 
in normal emission at a photon energy of 50 eV, and 
the indium Ad core level. The InP sample was in Ohmic 
contact with a Au metal surface cleaned by scratching, 
on which the position of the Fermi level was measured. 
With the STM we measured constant-current images at 
a variety of voltages. From these images the vacancy 
concentrations were determined as a function of the an- 
nealing time Jl3| as well as the spatial distribution of the 
occupied and empty density of states of the vacancies 
with atomic resolution. 

We first focus on the determination of the energy of 
the defect state. Here we concentrate on the charge trans- 
fer level of the surface vacancy, which is defined as the 
position of the surface Fermi level at which the defect 
changes its charge state. This energy level can be cal- 
culated by DFT and can be probed by measuring the 
band bending as shown below. We note that the energy 



of the charge transfer level might be very different from 
that of the ionization levels or the Kohn-Sham eigenval- 
ues within DFT. In particular, large deviations between 
these quantities might occur in systems with large atomic 
relaxations such as surface and bulk defects [ fl0| |. 

STM images reveal that during heat treatment of ini- 
tially defect-free cleavage surfaces the vacancy concentra- 
tion increases with increasing annealing time. Similarly, 
the position of the Fermi level, which on well-cleaved 
defect-free surfaces is close to the top of the valence 
band (VBM) for p-type samples, is found to shift toward 
midgap upon annealing, as measured by PES (see Fig. 
2). The filled diamonds (left axis) show the band bending 
at the surface, as determined through PES, whereas the 
open circles show the vacancy concentration (from STM, 
right axis) as a function of the annealing time. The band 
bending reaches a saturation value of 0.65 eV at vacancy 
concentrations of 5 x 10 12 cm~ 2 . Figure 2 also demon- 
strates that the Fermi level shift is directly correlated to 
the increase of the vacancy concentration. 

By combining the vacancy concentration and the band 
bending it induces, it is possible to determine the energy 
of the charge transfer level in the band gap E s d, because 
the charge per surface area in the surface layer Q ss is ex- 
actly compensated by the charge density per surface area 
in the depletion layer Q sc [l^j . The charge per surface 
area induced by a concentration n s d of +le charged P 
surface vacancies is [fLal 
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where Ep is the Fermi energy. The difference in energy of 
the charge transfer level and the Fermi level E sc i — Ep is 
given by (E sd - E sv ) - (E v - E sv ) - (E F - E v ) (for defini- 
tions, see inset in Fig. 2; E v and E sv are the positions of 
the valence band in the bulk and at the surface, respec- 
tively) where E v — E sv is the band bending eV s measured 
by photoelectron spectroscopy at the surface. The charge 
density Q sc in the depletion layer compensating Q ss is 



Qsc — \ 2eoe r ndo P kT 



-eV eV s . 
eXp( ^ )+ fcT- 1} 
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The energy difference between the Fermi energy and 
the valence band maximum in the bulk has been deter- 
mined for a nondegenerate InP crystal from the carrier 
concentration according to 



Ep- E v = -kT In 
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with m p being the mass of the holes. 

These equations permit a determination of the differ- 
ence in energy of the charge transfer level E sc i and the 
valence band maximum at the surface E sv using the va- 
cancy concentration dependence of the band bending (see 
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FIG. 2. The band bending (filled diamonds, left axis) and 

the vacancy concentration (empty circles right axis) on 

the InP(llO) surface as a function of the annealing time at 

160°C. The solid line is the vacancy concentration calculated 

on the basis of Eqs. (1) and (2) and a fit through the measured 

band bending values. The dashed line shows the saturation 

value. The right scale is chosen such that according to Eqs. 

(1) and (2) the vacancy concentration is shown on a scale 

linear proportional to the band bending. The inset shows the 

schematic drawing of the depletion region in case of a fully 

pinned surface. 



Fig. 2). We obtain an energy for the charge transfer level 
of the P vacancy of 0.75 ± 0.1 eV above the valence band 
maximum at room temperature. Note that the knowl- 
edge of the defect concentrations is crucial to determine 
the charge transfer level from a measurement of the band 
bending. 

At this stage we compare this energy value with re- 
sults from our total-energy density-functional theory cal- 
culations. Details of the calculational method used are 
described in Refs. |l7],[l8]]. Here we only increased the 
cutoff of the plane wave basis set to 15 Ry in the elec- 
tronic structure calculation to ensure convergence of the 
indium pseudopotential. We used a 2 x 4 and 6 layers 
thick supercell, which is large enough here, since the cal- 
culated bond lengths changed less than 0.01 A and the 
defect formation energy by less than 0.01 eV when in- 
creasing the size of the supercell. For the +le charged P 
vacancy the calculation yields two different equilibrium 
atomic configurations, whose total energies differ by only 
0.07 eV. The vacancy configuration with the lowest en- 
ergy exhibits a rebonding of one of the neighboring sur- 
face atoms with the indium in the second layer leading 
to a nonsymmetric atomic structure with respect to the 
(110) mirror plane of the defect-free (110) surface [see 
atomic structure superposed on Fig. 3(a)]. The calcu- 
lated indium- indium spacing of this configuration is 2.89 
A as compared to 2.98 A for the symmetric ground state 
structure of the neutral defect. A similar relaxation has 
also been calculated for the positively charged As va- 
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FIG. 3. Simulated STM images deduced from the calcu- 
lated local density of states 3 A above the surface of a +le 
charged P vacancy on InP(llO) surfaces for different atomic 
configurations (a) — (c) and experimental high resolution STM 
images of a P vacancy measured at (dl) —2.9 and (d2) +2.0 
V. The simulation shows the occupied (upper frames, energy 
range from VBM to VBM —1 eV) and empty (lower frames, 
energy range from VBM to VBM +2 eV) states for (al),(a2) a 
nonsymmetric rebonded; (bl),(b2) asymmetric nonrebonded; 
and (cl),(c2) the time average of the two mirror configuration 
of the nonsymmetric rebonded vacancy configuration. The 
respective atomic configurations are superposed on the simu- 
lated images. 



cancy on GaAs(llO) |2j,|l2|. From our calculations we 
find a charge transfer level (+/0) at 0.52 eV above the 
VBM for the nonsymmetric vacancy. For the energeti- 
cally less stable symmetric configuration this level is at 
0.45 eV. Both values deviate from the measured charge 
transfer level by roughly 0.2 eV. We attribute this devi- 
ation to the well-known underestimation of energy levels 
in local density approximation. An approximate upper 
limit of the error is the difference between the calculated 
and experimental band gap. Using our calculated value 
(1.26 eV) and the experimental one (1.42 eV) we get a 
systematic error of about 0.16 eV. We can therefore con- 
clude that the theoretical and experimental positions of 
the defect levels agree well within the error margins. We 
note, however, that the systematic error for the calcu- 
lated energies of the charge transfer levels is too large to 
identify the symmetry of the vacancy on the position of 
the defect level only. In the following we will show that 
an unambiguous identification can be obtained based on 
a comparison of the total energies and an analysis of the 
STM micrographs. 

We now focus on the atomic structure. STM images 
of +le charged P vacancies [Fig. 3(d)] always exhibit an 
apparently symmetric density of states with respect to 
the (110) mirror plane. The STM images thus suggest 
a symmetric vacancy structure although in our calcula- 
tions we find the nonsymmetric structure to be lower in 
energy. In order to clarify this apparent conflict, we have 
calculated in Fig. 3(a) — (c) simulations of occupied (up- 
per frames) and empty (lower frames) state STM images 
from the local density of states M] of different vacancy 
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configurations. Frames 3(al) and 3(a2) show the simu- 
lated STM images for the nonsymmetric case, whereas 
frames 3(bl) and 3(b2) show the same for the symmetric 
vacancy configuration. Both vacancy configurations lead 
to clearly distinguishable STM images: the nonsymmet- 
ric vacancy exhibits a pronounced asymmetry. The P 
atom next to the rebonded indium atom [upper atom in 
Fig. 3(a)] has a brighter dangling bond than the P atom 
at the opposite side of the vacancy and the rebonded in- 
dium atom has a weaker maximum in the empty states. 
Lateral displacements are also visible. These effects are 
purely electronic, because the calculated height differ- 
ence is only 0.11 and 0.07 A for the cations and anions, 
respectively. In contrast, the symmetric vacancy shows 
the expected symmetry with respect to the (110) mir- 
ror plane [Fig. 3(b)]. The images of the occupied states 
exhibit very little changes of the neighboring occupied 
dangling bonds [Fig. 3(bl)], and the separation of the 
two neighboring empty dangling bonds is considerably 
reduced [Fig. 3(b2)]. If we compare these simulations 
with the experimental STM images we find a symmetric 
density of states [Fig. 3(d)] in conflict with the energet- 
ically favorable nonsymmetric vacancy configuration in 
Fig. 3(a). However, also the symmetric vacancy configu- 
ration does not agree with the experimental STM images: 
The empty state STM images do not show any reduced 
separation of the two neighboring empty dangling bonds 
in contrast to the simulation of the symmetric vacancy 
structure [Fig. 3(b2)]. Thus none of the static equilib- 
rium vacancy configurations can be matched in detail to 
the experimental observation. 

This conflict can be resolved by closer inspection of the 
barrier between the two possible mirror configurations of 
the asymmetric vacancy. We estimated an upper limit 
for this barrier by mapping the energy along the direct 
reaction path between the symmetric and nonsymmet- 
ric vacancy configurations. On this basis we found an 
upper limit of 0.08 eV for the barrier between the two 
nonsymmetric configurations. The low value of the bar- 
rier implies that, at room temperature, the vacancy flips 
between the two nonsymmetric configurations at a rate 
significantly higher than the time resolution of the STM 
(1 to 10 Hz). We estimate a flip rate at room temperature 
of 0.1 THz using the transition state theory, the optical 
phonon frequency of 2 THz Jis|| as attempt frequency, 
and the calculated barrier. This concept of a thermally 
activated flip motion has been well established for dimers 
at the Si(001) surface pO|] . Based on these findings we 
interpret the STM image as a time average of the flipping 
vacancy. 

In order to test the model of a thermal flip motion we 
have performed STM image simulations based on this flip 
mechanism. In Fig. 3(c) we assumed that the vacancy 
flips such that the density of states probed by the STM is 
the average of the rebonded and nonrebonded sides of the 
vacancy (average of the two mirror configurations). The 



images agree very well with the experiment. In particu- 
lar, two major features — the nearly unchanged separa- 
tion of the neighboring empty dangling bonds [Fig. 3(c2)] 
and the depression of the neighboring occupied dangling 
bonds [Fig. 3(d)] — are well reproduced. Close inspec- 
tion of Figs. 3(c) and 3(d) clearly demonstrates that the 
images of the thermally flipping vacancy configuration 
explains the STM images best. Note that the somewhat 
higher intensity in the experimental empty state images 
is due to the local electrostatic potential arising from the 
positive chare Q and that the small intensity of empty 
density of states above the anions decays quickly into the 
vacuum such that the STM cannot probe it [|l|. 

In conclusion, we have demonstrated a methodology 
to determine the atomic structure and the energy of 
charge transfer levels in the band gap of surface defects, 
by combining scanning tunneling microscopy, photoclcc- 
tron spectroscopy, and density-functional theory calcula- 
tions. From linking STM and PES data we find a charge 
transfer level (+/0) 0.75 ±0.1 eV above the VBM for 
the example of P vacancies on InP(llO) in good agree- 
ment with DFT calculations. By comparing high reso- 
lution STM images with calculations, we determine the 
vacancy to exhibit a rebonded nonsymmetric structure. 
The vacancy flips thermally between its two mirror con- 
figurations. This localized vibrational mode leads to an 
apparently symmetric density of states in STM images, 
although the vacancy has a nonsymmetric atomic struc- 
ture. Similar measurements may help to understand on 
the atomic scale the properties of a wide range of mate- 
rials, where defects play a key role. 

We thank Th. Chasse for assistance in an initial ex- 
periment, W. Monch for valuable discussions, and K. H. 
Graf for technical support. This work was supported 
by the Bundesministerium fur Forschung und Technolo- 
gie under Grant No. 05SE8 OLA7 and by the Deutsche 
Forschungsgemeinschaft through Sonderforschungsbere- 
ich 296 TP A4 and A5 and Grant No. Ne428/2-l. 



[1] G. Lengel et ai, Phys. Rev. Lett. 72, 836 (1994). 
[2] S. B. Zhang and A. Zunger, Phys. Rev. Lett. 77, 119 
(1996). 

[3] H. Kim and J. R. Chelikowsky, Phys. Rev. Lett. 77, 1063 

(1996) ; Surf. Sci. 409, 435 (1998). 

[4] Ph. Ebert, K. Urban, and M. G. Lagally, Phys. Rev. Lett. 

72, 840 (1994). 
[5] J. Harper, G. Lengel, R. E. Allen, and M. Weimer, Phys. 

Rev. Lett. 79, 3312 (1997). 
[6] S. B. Zhang and A. Zunger, Phys. Rev. Lett. 79, 3313 

(1997) . 

[7] J. Harper, G. Lengel, R. E. Allen, and M. Weimer, Phys. 
Rev. Lett. 79, 3314 (1997). 



4 



[8] H. Kim and J. R. Chelikowsky, Phys. Rev. Lett. 79, 3315 
(1997). 

[9] S. Aloni, I. Nevo, and G. Haase, Phys. Rev. B 60, R2165 
(1999). 

[10] U. Scherz and M. Schemer, in Imperfections in III-V Ma- 
terials, edited by E. R. Weber (Academic Press, New 
York, 1993). 

[11] Ph. Ebert, Surf. Sci. Reports 33, 121 (1999). 

[12] The theoretical results for the positively charged arsenic 
vacancy on GaAs(llO) exhibit a very close agreement to 
those for the respective P vacancy on InP(llO). Thus 
both vacancies can be expected to show a thermal flip 
motion resulting in good agreement between STM simu- 
lation and experiment. 

[13] Ph. Ebert et al, Phys. Rev. B 51, 9696 (1995). 

[14] Ph. Ebert et al, Phys. Rev. B 53, 4580 (1996). 

[15] Ph. Ebert et al, Phys. Rev. Lett. 76, 2089 (1996). 

[16] W. Monch, Semiconductor Surfaces and Interfaces, 
Springer Series in Surface Sciences Vol. 26, edited by G. 
Ertl et al., (Springer, Berlin/Heidelberg, 1993). 

[17] G. Schwarz, A. Kley, J. Neugebauer, and M. SchefHer, 
Phys. Rev. B 58, 1392 (1998). 

[18] M. Bockstedte, A. Kley, J. Neugebauer, and M. Schef- 
fler, Comput. Phys. Commun. 107, 187 ( 1997); see also 
http:/ /www. fhi-berlin.mpg.de/th/fhimd/ 

[19] U. Schroder, J. Fritsch, and P. Pavone, Physica (Ams- 
terdam) 219B-220B, 434 (1996); H. Nienhaus and W. 
Monch, Surf. Sci. 328, L561 (195). 

[20] J. Dabrowski and M. Schemer, Appl. Surf. Sci. 56-58, 15 
(1992). 

[21] B. Engels et al, Phys. Rev. B 58, 7799 (1998). 



5 



